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The dynamical responses of Ising metamagnet (layered antiferromagnet) in 
the presence of a sinusoidally oscillating magnetic field are studied by Monte 
Carlo simulation. The time average staggered magnetization plays the role 
of dynamic order parameter. A dynamical phase transition was observed 
and a phase diagram was plotted in the plane formed by field amplitude and 
temperature. The dynamical phase boundary is observed to shrink inward as 
the relative antiferromagnetic strength decreases. The results are compared 
with that obtained from pure ferromagnetic system. The shape of dynamic 
phase boundary observed to be qualitatively similar to that obtained from 
previous meanfield calculations. 
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I. Introduction: 

The dynamic pliase transition[l], in pure Ising ferromagnet driven by 
oscillating magnetic field, became an interesting field of modern research 
in nonequilibrium statistical physics. The researchers paid much attention 
in last few years to study the dynamical phase transition in other mag- 
netic models also. The dynamic transition is studied [2] in kinetic spin-3/2 
Blume-Capel model. The multicritical behaviour was observed[3] in kinetic 
Blume-Emery-Griffith model. The dynamic transition was investigated in 
the classical Heisenberg model[4] with bilinear exchange anisotropy and in 
XY model [5]. The multiple dynamic phase transition was also observed [6] 
in classical anisotropic Heisenberg ferromagnet driven by polarised magnetic 
field. The existence of dynamic phase transition was found by few experi- 
mental studies in systems like ultrathin Co film on Cu(100)[7, 8] by surface 
magneto-optic Kerr effect. 

However all these studies mentioned above are done in simple ferromag- 
netic systems. Due to the presence of complicated exchange interactions, 
the dynamical responses of metamagnets driven by oscillating magnetic field 
may give rise to some interesting effects. Keeping this in mind, recently few 
researchers have taken interest to study the dynamical behaviours of meta- 
magnets driven by oscillating magnetic field. Few investigations are made in 
this front. Dynamic phase transition is studied in kinetic metamagnetic[9] 
spin-3/2 Blume-Capel model, in Ising metamagnets[10] etc. However, all 
these studies are mainly based on meanfield theory. 

In the meanfield calculations the spin fluctuations are ignored and the 
results do not have any informations of microscopic details. Moreover, the 
transition is studied only from the temperature variations of the order pa- 
rameters. The temperature variations of quantity like specific heat cannot 
be studied. One possible way to incorporate, the fluctuations as well as the 
temperature variations of specific heat, is to study this by Monte Carlo sim- 
ulations. As far as the auhor's knowledge is concerned, no such attempt has 
been made so far to study the dynamic phase transition, even in simple Ising 
metamagnet driven by oscillating magnetic field, by Monte Carlo simulation. 
Being motivated by these facts, the dynamic response of Ising metamagnet 
driven by oscillating magnetic field, is studied by Monte Carlo simulation, in 
this article. 

The paper is organised in the following manner. The Ising metamagnetic 
model and the Monte Carlo simulation scheme are discussed in the next 
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section (section-II). The numerical results are given in section-Ill. The paper 
ends with concluding remarks and summary in section-IV. 

II. Model and Simulation: 

The time dependent Hamiltonian (or energy) of a three dimensional Ising 
metamagnet (layered antiferromagnet) is represented as: 

H{t) = -Jp^pSiSj - JA^ASiSj - h{t)T,Si. (1) 

First term represents the in-plane ferromagnetic {Jp > 0) nearest neighbour 
spin-spin interaction energy. Second term provides the interaction energy 
coming from the antiferromagnetic {Ja < 0) interaction between two adjus- 
cent layers (since the nearest neighbour interactions arc considered only). 
Third term gives the spin (sj = ±l)-field {h{t)) interaction energy. Here the 
time dependent field is taken in the sinusoidal form, i.e., h{t) — hocos{ujt). 
The boundary condition is taken periodic in all directions. 

A cubic lattice of linear size L = 20 is considered. The dynamical 
evolution of the Ising spins are studied by Monte Carlo simulation using 
Metropolis[ll] single spin- flip scheme. The initial state of spin configura- 
tion is prepared by taking 50 percent of total number of spins (randomly 
selected) up (sj = -1-1). This corresponds to a high temperature paramag- 
netic phase. According to Metropolis [11] single spin flip dynamics, a spin, 
selected randomly, will flip (sj — )■ — Sj) with probability, 

W{si ^ -Si) = Min[l, bM-^H/KbT)] (2) 

where, AH is the change in energy due to spin flip, Kb is Boltzmann con- 
stant and T is the temperature. number of such random updates of spins 
is one Monte Carlo Step (MCS) and defines the unit of time in the present 
study. After a long time a dynamical steady state is achieved. The satisfac- 
tory steady values of the time averages of the dynamical physical quantities 
(defined in next section) ensures the achivement of steady state. All dy- 
namical quantities are calculated in the steady state for a fixed set of values 
of the temperature (T) and amplitude (ho) of the oscillating magnetic field. 
Then the temperature is reduced (by a small step) keeping the values of other 
parameters unchanged and a similar process (described above) is repeated. 
Here, the last spin configuration is used as the initial state for the present 
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value of temperature. In this way the temperature variations of all dynam- 
ical quantities are studied. Here, the temperature and field amplitudes are 
measured in the unit where, Kb — 1 and Jp — 1. 

III. Numerical results: 

The zero temperature configuration of this system is all spins are par- 
allel in all layeres and the adjuscent layers contains antiparallel spins. The 
system can be decomposed in two different sublattices (say A and B). So, 
alternate layeres form a sublattice. In a particular sublattice (A), the in- 
stantaneous magnetization is Mj^{t) = 2(SjSj)/L^. The time average (over a 
full cycle of the oscillating magnetic field) magnetization (for sublattice A) 
is Qa — ^ § MA{t)dt. In the present study the frequency of the oscillating 
magnetic field is taken u = 2%f = 27r x 0.01. The frequency (/) of the 
oscillating magnetic field is kept constant (/ = 0.01) throughout the study. 
So, one complete oscillation would require 100 MCS. Initially, data for 1200 
such cycles are discarded and average value of Qa is calculated from next 
300 cycles. It is checked that 300 number of cycles are sufficient to achieve 
the dynamically stable values of quantities. The dynamic order parameter 
Q B for other sublattice (B) is calculated in the same way. The staggered dy- 
namic order parameter Qs is calculated as the time average (over a full cycle 
of the oscillating magnetic field) of instantaneous staggered magnetization 
^ {MA-MB) y gg^g^ Qs = ^ §{MA{t) - MB{t))dt. The dynamical average en- 
ergy is also calculated as E — ^ f H{t)dt. The dynamical specific heat (C) 
is defined as C = ^ [12]. The temperature variations (with step AT = 0.05) 
of Qa, Qb, Qs, E and C are studied here, considering the amplitude (ho) 
and frequency (/) of the time dependent magnetic field as parameters in each 
case. 

Figure- 1 shows the temperature variations of dynamic quantities for fixed 
values of Jp = 1-0, Ja = —1.0, / = 0.01 and for two different values {ho = 2.5 
and Hq = 3.0). Figure-l(a) shows the temperature variations of the sublat- 
tice dynamic order parameters {Qa and Qb)- This variations clearly shows 
the two dynamic transitions near — 2.2 for ho — 2.5 and — 1.35 for 
ho — 3.0. One can easily visualize the temperature variations (not shown) of 
dynamic staggered order parameters {Qs)- These transitions becomes more 
pronounced from the plot of temperature variations of the derivatives (^^) 
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(in Figure- 1(b)). The sharp minima (eventually divergences in L — )■ oo limit) 
indicate the dynamic transition temperatures = 2.20 and T^, = 1.35 for 
ho — 2.5 and ho — 3.0 respectively. These two dynamic transition were recon- 
firmed and re-estimated independently from the studies of the temperature 
variations of dynamic specific heat (C)[12]. It may be mentioned here that 
the temperature variation of the dynamic specific heat is not studied by pre- 
vious meanfield calculations[10] These variations are shown in Figure-l(c). 
Here also, the dynamic specific heat (C) shows peaks (eventually divergences 
in the limit L — > oc) near — 2.20 and T^, — 1.35 for ho = 2.5 and ho = 3.0 
respectively. In this way, the transition temperatures for the dynamic phase 
transitions arc estimated over the range of values of the amplitudes of the 
oscillating magnetic field. Thus for fixed values of Ja = —1.0, Jp = 1-0 and 
/ = 0.01, the dynamic phase boundary is obtained. 

To study the variations of dynamic phase boundary with the relative 
strength of the antiferromagnetic interactions {Ja), the similar investigations 
are made for Ja = —0.5 and Jp = 1.0. In this case, the relative antiferro- 
magnetic strength is reduced. Figure-2(a) shows the temperature variations 
of the derivative of staggered dynamic order parameters, i.e., Like the 
previous cases, the shows sharp minima indicating the dynamic phase 
transitions near Td = 1.20 for ho — 2.0 and near Td = 0.45 for ho = 3.0. The 
same transition temperatures were recalculted from the temperature varia- 
tions of dynmaic specific heat (C). The dynamic specific heat C shows sharp 
maxima near the same transition points indicating the dynamic transitions. 
These results are shown in Figure-2(b). In this way, the entire dynamic phase 
boundary is obtained. 

For a comparison, the dynamic transitions are studied for a pure ferro- 
magnetic system. The transition temperatures = 1.6 and = 0.8 are 
obtained from the temperature variation of the derivative of dynamic order 
parameter ^, for Q = ^ § M{t)dt, where M{t) = [MA(t) + Ms(t)] is the to- 
tal instantaneous magnetization obtained by piitting Jp — 1.0 and J a = 1.0 
in equation (1). The dynamic transition temperatures are estimated from 
the sharp minima (shown in Figure-3) of The results obtained here, are 
Td — 1.6 for ho — 2.0 and — 0.8 for ho = 3.0. Hence the data for entire 
dynamic phase bounadry are obtained. 

The comprehensive results of the dynamic phase transitions are shown 
as the dynamic phase boundary in Figure-4. From the phase diagrams it is 
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clear that the shapes of dynamic phase boundaries for metamagnetic dynamic 
phase transitions are distinctly different from that for pure ferromagnets. 
Additionally, the dynamic phase boundary shrinks inward as the relative 
strength of antiferromagnetic interaction decreases. These phase diagrams 
for metamagnetic dynamic phase transitions agrees qualitatively well with 
those obtained from meanfield calculations [10]. 

IV. Concluding remarks: 

The dynamic phase transition in Ising metamagnet driven by oscillating 
magnetic field is studied by Monte Carlo simulation. As far as the author's 
knowledge is concerned, this is the first Monte Carlo study of dynamic phase 
transition for this type of metamagnetic model. This study differes from the 
previous meanfield study [10] in the following aspects: Firstly, the dynamic 
phase transition was reconfirmed here from the temperature variation of dy- 
namic specific heat [12] , which was not done in meanfield study [10]. Secondly, 
in meanfield study, the multicritical behaviours are observed and the tricrit- 
ical point was located on the phase boundary at lower temperatures. Here, 
in the Monte Carlo study, the probability of spin fiip was calculated from the 
Gibb's distribution which provides good results in high temperature region. 
In the low temperature, region the quantum fiuctuation will be excited and 
Quantum Monte Carlo will be the better method of study. For this reason, 
the phase boundary for very low temperature is not drawn. 

In the present study, the order of transition is not mentioned. As a result 
one cannot be able to guess about the multicritical behaviour. Here, only the 
dynamic transition is observed. It may be mentioned here, that the existence 
of tricritical point was found[13] on the dynamic phase boundary by Monte 
Carlo simulation in pure Ising ferromagnet. However, this was disproved[14] 
later from the Monte Carlo study with larger lattice size and improved statis- 
tics. So, it is better not to mention about multicritical behaviour from the 
present Monte Carlo study of such small system size {L = 20). 

The shapes of the dynamic phase boundary depend on the relative anti- 
ferromagnetic strength as well the the frequencies of the oscillating magnetic 
field. In the present case, this variation with frequency is not studied. It 
may be an interesting study. 

Above all, this dynamic phase transition in Ising metamagnets gives in- 
teresting result, theoretically. The physical reason of the dynamic phase 
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transition in metamagnet may be stated qualitatively as follows: for small 
values of the amplitude of the field and at low temperature, the instanta- 
neous staggered magnetization oscillates about a nonzero value, giving rise 
to ordered phase. However, for large value of field amplitude and at high 
temperature, the instantaneous staggered magnetization oscillates symmet- 
rically about a value very close to zero, leading to disordered phase. One has 
to see the hysteretic loop (m — h) [1] to get the clear idea about the physical 
interpretation of the dynamic phase transition. The asymmetric m — h loop 
gives the dynamically ordered phase and the symmetric m — h loop gives 
disordered phase. In the case of metamagnet, this m should be the staggered 
magnetization. This theoretical investigation is an apeal to the experimen- 
talists to study the dynamic phase transition in anisotropic metamagnet, 
like FeBr2. Experimentally the existence of dynamic transition was found 
[7, 8] in Co film (ferromagnetic) on Cu surface (at room temperature) by 
surface magneto optic Kerr effect. Let us hope that the experimental study 
of dynamic phase transition in metamagnetic systems like FeBr2 driven by 
oscillating magnetic field, will bring some new physics in future. 
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Fig-1. The temperature variations of (a) sublattice dynamic order parame- 
ters Qa ((•) for ho = 3.0 and (o) for ho = 2.5)and Qb ((A) for ho = 3.0 and 
(o) for ho = 2.5), (b) the derivative (^?) of staggered dynamic order param- 
eter (Qs) ((o) for ho — 3.0 and (o) for ho — 2.5) and (c) the dynamic specific 
heat (C) ((o) for ho — 3.0 and (o) for ho — 2.5). The continuous hues in all 
cases joining the data points act as guides to the eye. Here, — —1 x Jp. 
The Td{ho = 2.5) = 2.2 and Td{ho = 3.0) = 1.35. 
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Fig-2. The temperature variations of (a) for = 3.0 (o) and = 2.0(o) 
and (b) C for = 3.0 (o) and ho = 2.0 (o). The continuous hues in all 
cases joining the data points act as guides to the eye. Here, Ja = —0.5Jf- 
Td{ho = 3.0) = 0.45 and Td{ho = 2.0) = 1.20. 
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Fig-4. The phase diagram of dynamic phase transitions. Metamagnetic 
dynamic phase bounadries (i) (o) for Ja = —0.5 and Jp = 1.0 and (ii) (o) for 
Ja — —1-0 and Jp — 1.0. The (•) represents the dynamic phase bounadry 
of pure ising ferromagnets (J^ = Ji;' = 1.0). 
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